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Abstract 

The relativistic equation of state (EOS) of the Van-der-Waals gas is suggested and 
analyzed. In contrast to the usual case, the Lorentz contraction of the sphere's volume 
is taken into account. It is proven that the suggested EOS obeys the causality in the 
limit of high densities, i.e., the value of sound velocity of such a media is sublumi- 
nar. The pressure obtained for the high values of chemical potential has an interesting 
kinetic interpretation. The suggested EOS shows that for high densities the most prob- 
able configuration corresponds to the smallest value of the relativistic excluded volume. 
In other words, for high densities the configurations with the collinear velocities of the 
neighboring hard core particles are the most probable ones. This, perhaps, may shed 
light on the coalescence process of any relativistic hard core constituents. 



Key words: Equation of state, hard spheres, relativistic Van-der-Waals model 



1 



1. Introduction 



The van der Waals (VdW) excluded volume model is successfully used to describe the 
hadron yields measured in relativistic nucleus-nucleus collisions (see e.g. [1, 2] and refer- 
ences therein). This model treats the hadrons as hard core spheres and, therefore, takes 
into account the hadron repulsion at short distances. In a relativistic situation one should, 
however, include the Lorentz contraction of the hard core hadrons. Recently, both the con- 
ventional cluster and the virial expansions were generalized to the momentum dependent 
inter-particle potentials, accounting for the Lorentz contracted hard core repulsion [3] and 
the derived equation of state (EOS) was applied to describe hadron yields observed in rel- 
ativistic nuclear collisions [4]. The VdW equation obtained in the traditional way leads to 
the reduction of the second virial coefficient (analog of the excluded volume) compared to 
nonrelativistic case. However, in the high pressure limit the second virial coefficient remains 
finite. This fact immediately leads to the problem with causality in relativistic mechanics - 
the speed of sound exceeds the speed of light [5] . 

The influence of relativistic effects on the hard core repulsion may be important for a 
variety of effective models of hadrons and hadronic matter such as the modified Walecka 
model [6], various extensions of the Nambu-Jona-Lasinio model [7], the quark- meson cou- 
pling model [8], the chiral SU(3) model [9] e.t.c. Clearly, the relativistic hard core repulsion 
should be important for any effective model in which the strongly interacting particles have 
the reduced values of masses compared to their vacuum values because with lighter masses 
the large portion of particles becomes relativistic. Nevertheless, the relativistic hard core 
repulsion was, so far, not incorporated into these models due to the absence of the required 
formalism. 

The Lorentz contraction of rigid spheres representing the hadrons may also be essential 
at high particle densities which can be achieved at modern colliders. Very recently it was 
understood that in the baryonless deconfined phase above the cross-over temperature T c some 
hadrons may survive up to large temperatures like 3T C [10, 11, 12], and that above T c there 
may exist bound states [13] and resonances [14]. Moreover, an exactly solvable statistical 
model of quark-gluon bags with surface tension [15] indicates that above the cross-over 
transition [12] the coexistence of hadronic resonances with QGP may, in principle, survive 
up to infinite temperature. Thus, above T c the relativistic effects of the hard core repulsion 
can be important for many hardonic resonances and hadron-like bound states of quarks, 
especially, if their masses are reduced due to chiral symmetry restoration. 

Also the VdW EOS, which obeys the causality condition in the limit of high density and 
simultaneously reproduces the correct low density behavior, adds a significant theoretical 
value because such an EOS had not yet been formulated during more than a century of the 
special relativity. This work is devoted to the investigation of the necessary assumptions to 
formulate such an equation of state. 

The work is organized as follows. In Sect. 2 a summary of both the cluster and virial 
expansion for the Lorentz contracted rigid spheres is given. It is shown that the VdW ex- 
trapolation in relativistic case is not a unique procedure. Therefore, an alternative derivation 
of the VdW EOS is considered there. The high pressure limit is studied in details in Sect. 
3. It is shown that the suggested relativistic generalization of the earlier approach [3] obeys 
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the causality condition. The conclusions are given in the last section. 



2. Relativization of the van der Waals EOS 



The excluded volume effect accounts for the blocked volume of two spheres when they 
touch each other. If hard sphere particles move with relativistic velocities it is necessary to 
include their Lorentz contraction in the rest frame of the medium. The model suggested 
in Ref. [16] is not satisfactory: the second virial coefficient a 2 = 4-y of the VdW excluded 
volume model is confused there with the proper volume v Q of an individual particle - the 
contraction effect is introduced for the proper volume of each particle. In order to get the 
correct result it is necessary to account for the excluded volume of two Lorentz contracted 
spheres. 

Let Yi and ij be the coordinates of the i-th and j-th Boltzmann particle, respectively, 
and kj and kj be their momenta, r„ denotes the unit vector = r^/lr^l (r^ = |r» — r j- j ) . 
Then for a given set of vectors (f^, kj, kj) for the Lorentz contracted rigid spheres of radius 
R Q there exists the minimum distance between their centers ^(f y-; kj, kj) = min|i"jj|. The 
dependence of the potentials Uij on the coordinates r,, ij and momenta kj,kj can be given 
in terms of the minimal distance as follows 



kj; Tj, kj) < 



f o, 

< 




Jij] kj, kj) 






Jij', kj, kj) 



(1) 



The general approach to the cluster and virial expansions [17] is valid for this momentum 
dependent potential, and in the grand canonical ensemble it leads to the transcendental 
equation for pressure [3] 

p(T,p) = Tp t (T) exp ^ 2P ^j = p id (T, /i - a 2 p) , (2) 

with the second virial coefficient 

a 2 {T) = -I J ^ 6 e t u(ki,k 2 ), (3) 

^(k!,k 2 ) = l - J dT 12 Q(ru(T&MM) ~ |ri 2 |) , (4) 

E(k) 

where the thermal density is defined as follows p t (T) = g J j0js^ T , degeneracy as g, and 
f(ki,k2) denotes the relativistic analog of the usual excluded volume for the two spheres 
moving with the momenta ki and k 2 and, hence, the factor 1/2 in front of the volume integral 
in (4) accounts for the fact that the excluded volume of two moving spheres is taken per 
particle. 

In what follows we do not include the antiparticles into consideration to keep it simple, 
but this can be done easily. Then the pressure (2) generates the following particle density 

n(T,») = = f MT) V, (5) 

^ 1 + er p t (T)a 2 (T) T (l + eT p t (T)a 2 (T) 



3 



which in the limit of high pressure p(T, n) — > oo gives a limiting value of particle density 
n(T,n) - a^(T). 

A form of Eq. (2) with constant a<i was obtained for the first time in Ref. [6]. The new 
feature of Eq. (2) is the temperature dependence of the excluded volume a 2 (T) (3) which 
is due to the Lorentz contraction of the rigid spheres. This is a necessary and important 
modification which accounts for the relativistic properties of the interaction. It leads, for 
instance, to a 50 % reduction of the excluded volume of pions already at temperatures 
T = 140 MeV [3]. 

The calculation of the cluster integral in relativistic case is more complicated because 
each sphere becomes an ellipsoid due to the Lorentz contraction and because the relativistic 
excluded volume strongly depends not only on the contraction of the spheres, but also on 
the angle between the particle 3-velocities. Therefore, in Appendix A we give a derivation 
of a rather simple formula for the coordinate space integration in a 2 which is found to be 
valid with an accuracy of a few percents for all temperatures. Its simplicity enables us to 
perform the angular integrations in a 2 (T) analytically and obtain 

irp\ 0lv o ( 74p s \ r dk m _| 

The expression for the coefficient a 1/1.065 is given in Appendix A by Eq. (51). Using 
this result it is easy to show that in the limit of high temperature T ^> m the ratio of 
the scalar density p s (T) to the thermal density pt{T) in (6) vanishes and the second virial 
coefficient approaches the constant value: 



which is about || times smaller compared to the value of the nonrelativistic excluded volume, 
and, hence, is surprisingly very close to the dense packing limit of the nonrelativistic hard 
spheres. Similarly to the nonrelativistic VdW case [5] this leads to the problem with causality 
at very high pressures. Of course, in this formulation the superluminar speed of sound 
should appear at very high temperatures which are unreachable in hadronic phase. Thus the 
simple "relativization" of the virial expansion is much more realistic than the nonrelativistic 
description used in Refs. [1, 2], but it does not solve the problem completely. 

The reason why the simplest generalization (2) fails is rather trivial. Eq. (2) does 
not take into account the fact that at high densities the particles disturb the motion of 
their neighbors. The latter leads to the more compact configurations than predicted by 
Eqs. (2-4), i.e., the motion of neighboring particles becomes correlated due to a simple 
geometrical reason. In other words, since the iV-particle distribution is a monotonically 
decreasing function of the excluded volume, the most probable state should correspond to 
the configurations of smallest excluded volume of all neighboring particles. This subject is, 
of course, far beyond the present paper. Although we will touch this subject slightly while 
discussing the limit /i/T >> 1 in Sect. 3, our primary task here will be to give a relativistic 
generalization of the VdW EOS, which at low pressures behaves in accordance with the 
relativistic virial expansion presented above, and at the same time is free of the causality 
paradox at high pressures. 
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In our treatment, we will completely neglect the angular rotations of the Lorentz con- 
tracted spheres because their correct analysis can be done only within the framework of 
quantum scattering theory which is beyond the current scope. However, it is clear that the 
rotational effects can be safely neglected at low densities because there are not so many 
collisions in the system. At the same time the rotations of the Lorentz contracted spheres 
at very high pressures, which are of the principal interest, can be neglected too, because 
at so high densities the particles should be so close to each other, that they must prevent 
the rotations of neighboring particles. Thus, for these two limits we can safely ignore the 
rotational effects and proceed further on like for the usual VdW EOS. 

Eq. (2) is only one of many possible VdW extrapolations to high density. As in non- 
relativistic case, one can write many expressions which will give the first two terms of the 
full virial expansion exactly, and the difference will appear in the third virial coefficient. In 
relativistic case there is an additional ambiguity: it is possible to perform the momentum 
integration, first, and make the VdW extrapolation next, or vice versa. The result will, 
evidently, depend on the order of operation. 

As an example let us give a brief "derivation" of Eq. (2), and its counterpart in the grand 
canonical ensemble. The two first terms of the standard cluster expansion read as [17, 3] 

p = T Pt (T) er (l-a 2Pt {T)eT^ . (8) 

i£ 

Now we approximate the last term on the right hand side as pt{T)eT m Then we 
extrapolate it to high pressures by moving this term into the exponential function as 

p « T Pt (T) er (l - a 2 |) « T p t (T) exp (^=^) . (9) 

The resulting expression coincides with Eq. (2), but the above manipulations make it simple 
and transparent. Now we will repeat all the above steps while keeping both momentum 
integrations fixed 

Tg 2 eT r dk! dk 2 v(k l ,k 2 )p\ 
P ~ — 7^r / , n vi , n ^ e T 1 ~ 



Pt (T) J (27r)3(27r)3 V T J 

To 2 r dki dk 2 ^-v(^iM)p-E(k 1 )-E(k 2 ) 

1 e t . (10) 



Pt (T)J (27T)3(27T) 

The last expression contains the relativistic excluded volume (4) explicitly and, as can 
be shown, is free of the causality paradox. This is so because at high pressures the main 
contribution to the momentum integrals corresponds to the smallest values of the excluded 
volume (4). It is clear that such values are reached when the both spheres are ultrarelativistic 
and their velocities are collinear. 

With the help of the following notations for the averages 

g* r dk, dk 2 vtoMp + EM + Ete) 

{{ " ~ PKT)J (2vr)3 (2vr)3 ° ' (12) 
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we can define all other thermodynamic functions as 



n ( T 'N = a. = — 7 5 V' ( 13 ) 



dfx r(l + efp t (T)(Mk 1 ,k 2 )))) 



9T T T l + eT Pt (T)((^(k 1 ,k 2 ))) 

2ef Pf (T)((E))-[/i+(E)]pT- 1 
e(T,/x) = Ts(T,fi) + /xn(T,/x) -p(T,/x) = s . (15) 

l + eTp t (T)(( v (k 1 ,k 2 ))) 

Here n(T, /x) is the particle density, while s(T, /x) and e(T, /x) denote the entropy and energy 
density, respectively. 

In the low pressure limit Apv Q T~ l <C 1 the corresponding exponent in (12) can be 
expanded and the mean value of the relativistic excluded volume can be related to the 
second virial coefficient a 2 (T) as follows 

(Mki,k 2 ))>^2(T) - ^(kx.ka))), (16) 

which shows that at low pressures the average value of the relativistic excluded volume 
should match the second virial coefficient a 2 (T), but should be smaller than a 2 (T) at higher 
pressures and this behavior is clearly seen in Fig. 1. 

A comparison of the particle densities (5) and (13) shows that despite the different 
formulae for pressure the particle densities of these models have a very similar expression, 
but in (13) the second virial coefficient is replaced by the averaged value of the relativistic 
excluded volume ((v (k 1; k 2 ))). Such a complicated dependence of the particle density (13) 
on T and /x requires a nontrivial analysis for the limit of high pressures. 

To analyze the high pressure limit p — > oo analytically we need an analytic expression for 
the excluded volume. For this purpose we will use the ultrarelativistic expression derived in 
the Appendix A: 

n i \ v¥ 2 rel (R,R) v Q ( m m \ / 2 fQ v \\ 2 3x; . _ , 

„ (kl , k2 ) K ^^-^_ + _)(i W ( T )) + _ m( e.). (i7) 

As usual, the total excluded volume vY{ el {R,R) is taken per particle. Eq. (17) is valid for 
< & v < |; to use it for | < Q v < tt one has to make a replacement Q v — > tt — Q v 
in (17). Here the coordinate system is chosen in such a way that the angle Q v between 
the 3-vectors of particles' momenta ki and k 2 coincides with the usual spherical angle of 
spherical coordinates (see Appendix A). To be specific, the OZ-axis of the momentum space 
coordinates of the second particle is chosen to coincide with the 3-vector of the momentum 
ki of the first particle. 

The Lorentz frame is chosen to be the rest frame of the whole system because otherwise 
the expression for pressure becomes cumbersome. Here v Q stands for the eigen volume of 
particles which, for simplicity, are assumed to have the same hard core radius and the same 
mass. 
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Fig. 1. Comparison of the exact value of the second virial coefficient a,2(T)/ 0,2(0) (solid 
curve) with the averaged value of the relativistic excluded volume a((v (ki, k 2 )))/(a 2 (0)) 
(dashed curve) given by Eq. (17) fot fi = 0. The normalization coefficient a ~ 1/1.065 (51) 
is introduced to reproduce the low density results. 

Despite the fact that this equation was obtained for ultrarelativistic limit, it is to a within 
few per cent accurate in the whole range of parameters (see Fig. 1 and Appendix A for the 
details), and, in addition, it is sufficiently simple to allow the analytical treatment. 



3. High Pressure Limit 

As seen from the expression for the relativistic excluded volume (17) , for very high 
pressures only the smallest values of the relativistic excluded volume will give a non- vanishing 
contribution to the angular integrals of thermodynamic functions. This means that only G„- 
values around and around 7r will contribute into the thermodynamic functions (see Fig. 
2). Using the variable x = sin 2 (B„/2), one can rewrite the k 2 angular integration as follows 



Ie(h) = / * l"dx e -(--('-i) a +*V^) 



(18) 



with A = ^; B= 3 -A; C = + ^) , (19) 

where we have accounted for the fact that the integration over the polar angle gives a factor 
2tt and that one should double the integral value in order to integrate over a half of the Q v 
range. 
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Fig. 2. Comparison of the relativistic excluded volumes for highly contracted spheres. In 
the left panel the loner dashed curve corresponds to MM = 2 and E ^ = 10 whereas the 

* ° 1 m m 

short dashed curve is found for MM = 5 anc i MM _ ^he corresponding values in the 

mm 1 

right panel are MM = 10 MM = 10 (long dashed curve) and MM = 10 MM = 10 o 
(short dashed curve). It shows that the excluded volume for 9^ close to | is finite always, 
while for the collinear velocities the excluded volume approaches zero, if both spheres are 
ultrarelativistic. 



Since C < 2 in (19) is a decreasing function of the momenta, then in the limit A ^> 1 one 
can account only for the sjx dependence in the exponential in (18) because it is the leading 
one. Then integrating by parts one obtains 



Ie(ki 



4 



d k 2 k2 _ AC 



(27T) 



0.5 



dx e 



(20) 



Applying the above result to the pressure (10), in the limit under consideration one finds 
that the momentum integrals are decoupled and one gets the following equation for pressure 



P{T, 



16 T 3 



9v%p 2 p t (T) 
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dkk 2 



E(k) 



2v m 
TE(k) 



-1 2 



V 



[21) 



Now it is clearly seen that at high pressures the momentum distribution function in (21) may 
essentially differ from the Boltzmann one. To demonstrate this we can calculate an effective 
temperature by differentiating the exponential under the integral in (21) with respect to 
particle's energy E: 



Teff(E) 



_d_ 

dE 



E 
T 



2v m 
TE 



P 



T 



1 - 



2_Uo_m 
E 2 P 



(22) 



Eq. (22) shows that the effective temperature T e ff(E — > 00) = T may be essentially lower 
than that one at E = m. In fact, at very high pressures the effective temperature T e ff(m) 
may become negative. 
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A sizable difference between T e ff values at high and low particle energies mimics the 
collective motion of particles since a similar behavior is typical for the transverse energy 
spectra of particles having the collective transverse velocity which monotonically grows with 
the transverse radius [18, 19]. However, in contrast to the true collective motion case [18, 19] , 
the low energy T e ff (22) gets higher for smaller masses of particles. Perhaps, such a different 
behavior of low energy effective temperatures can be helpful for distinguishing the high 
pressure case from the collective motion of particles. 

Our next step is to perform the gaussian integration in Eq. (21). Analyzing the function 

F = 2\nk-^-A-^- (23) 
T E(k) v ; 

for A 1, one can safely use the ultrarelativistic approximation for particle momenta 
k ps E(k) — > oo. Then it is easy to see that the function F in (23) has an extremum at 



dF 2 1 A m „ Am m Am 

oE = E-f + A w> = =* E = E ~ jT7srz-r T (f + — + 1 ' ' (24) 

which turns out to be a maximum, since the second derivative of F (23) is negative 

2 2A^<0. (25) 



dE 2 



e=e* {E*y (E*y 



There are two independent ways to increase pressure: one can increase the value of 
chemical potential while keeping temperature fixed and vice versa. We will consider the 
high chemical potential limit /i/T>l for finite T first, since this case is rather unusual. In 
this limit the above expressions can be simplified further on 



E* ps ^2 



mv p. 



d 2 F 



dE 2 



e=e* T y / 2mv Q p 



(26) 



Here in the last step we explicitly substituted the expression for A. Performing the gaussian 
integration for momenta in (21), one arrives at 

rdkk 2 2v m (E*) 2 , 2E* 

fWf* T -^f^v^e— (27) 
which leads to the following equation for the most probable energy of particle E* 

E'^DT'e^, (28) 

9ir 6 p t (T) 

As one can see, Eq. (28) defines pressure of the system. Close inspection shows that the 
high pressure limit can be achieved, if the exponential in (28) diverges much slower than 
fi/T. The latter defines the EOS in the leading order as 

2 

E* ps H , =}► p Pd — d . (29) 

4' y 32mv K ' 
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The left hand side equation above demonstrates that in the fj,/T >> 1 limit the natural 
energy scale is given by a chemical potential. This is a new and important feature of the 
relativistic VdW EOS compared to the previous findings. 

The right hand side Eq. (29) allows one to find all other thermodynamic functions in 
this limit from thermodynamic identities: 

2p 

s « 0, — , e = Ts + fin — p m p . (30) 

[i 

Thus, we showed that for fx/T >> 1 and finite T the speed of sound c s in the leading order 
does not exceed the speed of light since 



2 d P 

C ' = Te 



= T C = 1 - (31) 

s/n as 



From Eq. (28) it can be shown that the last result holds in all orders. 

It is interesting that the left hand side Eq. (26) has a simple kinetic interpretation. 
Indeed, recalling that the pressure is the change of momentum during the collission time one 
can write (24) as follows (with E* = k*) 

= i^f = 2^ 1 (32) 

P 2mv tiR 2 8R 2' 1 ' 

In the last result the change of momentum during the collision with the wall is 2 A;*, which 
takes the time 3 8 ^°, . The latter is twice of the Lorentz contracted height (4/3i? ) of the 
cylinder of the base ttRI which is passed with the speed v*. Here the particle velocity v* and 
the corresponding gamma-factor 7* are defined as v *7* = k*/m. The rightmost factor 1/2 in 
(32) accounts for the fact that only a half of particles moving perpendicular to the wall has 
the momentum —k*. Thus, Eq. (32) shows that in the limit under consideration the pressure 
is generated by the particle momenta which are perpendicular to the wall. This, of course, 
does not mean that all particles in the system have the momenta which are perpendicular to 
a single wall. No, this means that in those places near the wall where the particles' momenta 
are not perpendicular (but are parallel) to it, the change of momentum 2k* is transferred 
to the wall by the particles located in the inner regions of the system whose momenta are 
perpendicular to the wall. Also it is easy to deduce that such a situation is possible, if the 
system is divided into the rectangular cells or boxes inside which the particles are moving 
along the height of the box and their momenta are collinear, but they are perpendicular to 
the particles' momenta in all surrounding cells. Note that appearing of particles' cells is a 
typical feature of the treatment of high density limit [20] and can be related to a complicated 
phase structure of nuclear matter at very low temperatures [21]. 

Of course, inside of such a box each Lorentz contracted sphere would generate an excluded 
volume which is equal to a volume of a cylinder of height and base itR^. This cylinder, 
of course, differs from the cylinder involved in Eq. (32), but we note that exactly the 
hight 4% is used in the derivation of the ultrarelativistic limit for the relativistic excluded 
volume (50) (see Appendix A for details). Thus, it is very interesting that in contrast to 
nonrelativistic case the relativistic excluded volume ^p 2 - which enters into Eq. (32) is only 
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33 % smaller than the excluded volume of ultrarelativistic particle at high pressures. 
Also it is remarkable that the low density EOS extrapolated to very high values of the 
chemical potential, at which it is not supposed to be valid at all, gives a reasonable estimate 
for the pressure at high densities. 

Another interesting conclusion that follows from this limit is that for the relativistic VdW 
systems existing in the nonrect angular volumes the relativistic analog of the dense packing 
may be unstable. 

The analysis of the limit T//i > 1 and finite \x also starts from Eqs. (21)-(24). The 
function F from (23) again has the maximum at E* = E(k*) = k* defined by the right hand 
side Eq. (24). Now the second derivative of function F becomes 



d 2 F 



dE 2 



m 



2^1 + 4^ 



-2A-— = - * ' . (33) 



E=E , (E*) 2 (E*) 3 (E*) 



This result allows one to perform the gaussian integration for momenta in (21) for this limit 
and get 

_9 ( i _i_ Ahl\ 1 

rdkk 2 £W 2v m (E*) 3 e \ t J / Am\ 
where the auxiliary integral 1^ is defined as follows 

+oo 

J € (s) = j die-?. (35) 



i 



The expression (34) can be also used to find the thermal density pt{T) in the limit T — > oo 
by the substitution A = 0. Using (34), one can rewrite the equation for pressure (21) as the 
equation for the unknown variable z = Am/T = 2v ^J ip 

z 3 *eT<p(z), 4>{z) = —. ^ V^. (36) 

(Syre 2 -^- 1 ) I ( (l)(l + z)2 

Before continuing our analysis further on, it is necessary to make two comments con- 
cerning Eq. (36). First, rewriting the left hand side Eq. (36) in terms of pressure, one can 
see that the value of chemical potential is formally reduced exactly in three times. In other 
words, it looks like that in the limit of high temperature and finite \x the pressure of the 
relativistic VdW gas is created by the particles with the charge being equal to the one third 
of their original charge. Second, due to the nonmonotonic dependence of <f>(z) in the right 
hand side Eq. (36) it is possible that the left hand side Eq. (36) can have several solutions 
for some values of parameters. Leaving aside the discussion of this possibility, we will further 
consider only such a solution of (36) which corresponds to the largest value of the pressure 
(21). 

Since the function <f>(z) does not have any explicit dependence on T or //, one can establish 
a very convenient relation 

dz _n dz 
8T~ T dfj, 
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between the partial derivatives of z given by the left hand side Eq. (36). Using (37), one can 
calculate all the thermodynamic functions from the pressure p = f3 T 2 z (with e (3m v )~ r ) 
as follows: 

n « (3T 2 ^, (38) 



2Tz + T 2 ^ 



2p -J n , (39) 



e = Ts + fin — p^p. (40) 

The last result leads to the causality condition (31) for the limit T///> 1 and finite fx. 

In fact, the above result can be extended to any /i > — oo and any value of T satisfying 
the inequality 

E* « T(Vl + ^+ l) > m, (41) 
which is sufficient to derive Eq. (36). To show this, it is sufficient to see that for z = there 

iL iL 

holds the inequality z 3 < eT$(z\ which changes to the opposite inequality z 3 > eT${z) for 
z = oo. Consequently, for any value of \i and T satisfying (41) the left hand side Eq. (36) 
has at least one solution z* > for which one can establish Eqs. (37)-(40) and prove the 
validity of the causality condition (31). 

The model (10) along with the analysis of high pressure limit can be straightforwardly 
generalized to include several particle species. For the pressure p(T, {/^}) of the mixture of 
iV-species with masses rrii (i = {1, 2, .., N}), degeneracy c/j, hard core radius Ri and chemical 
potentials fa is defined as a solution of the following equation 

where the relativistic excluded volume per particle of species i (with the momentum ki) 
and j (with the momentum k 2 ) is denoted as i; ij -(k 1 ,k 2 ), E^ki) = \Jk\ + mf and Ej(k 2 ) = 



+ mj are the corresponding energies, and the total thermal density is given by the 
expression 

Ptot (T,M)= /jj^E t . (43) 



The excluded volume %(ki, k 2 ) can be accurately approximated by a vY{ el (Ri, Rj)/2 defined 
by Eqs. (50) and (51). 

The multicomponent generalization (42) is obtained in the same sequence of steps as the 
one-component expression (10). The only difference is in the definition of the total thermal 
density (43) which now includes the chemical potentials. Note also that the expression (42) 
by construction recovers the virial expansion up to the second order at low particle densities, 
but it cannot be reduced to any of two extrapolations which are suggested in [22] and [23] 
for the multicomponent mixtures and carefully analyzed in Ref. [4]. Thus, the expression 
(42) removes the non-uniqueness of the VdW extrapolations to high densities, if one requires 
a causal behavior in this limit. 
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4. Concluding Remarks 



In this work we proposed a relativistic analog of the VdW EOS which reproduces the 
virial expansion for the gas of the Lorentz contracted rigid spheres at low particle densities 
and is causal at high densities. As one can see from the expression for particle density 
(13) and from the corresponding relation for effective temperature (22) the one-particle 
momentum distribution function has a more complicated energy dependence than the usual 
Boltzmann distribution function, which would be interesting to check experimentally. Such 
a task involves considerable technical difficulties since the particle spectra measured in high 
energy nuclear collisions involve a strong collective flow which can easily hide or smear 
the additional energy dependence. However, it is possible that such a complicated energy 
dependence of the momentum spectra and excluded volumes of lightest hardons, i.e. pions 
and kaons, can be verified for highly accurate measurements, if the collective flow is correctly 
taken into account. The latter adjustment is tremendously complex because it is related to 
the freeze-out problem in relativistic hydrodynamics [24] or hydro-cascade approach [25]. 
Another possibility to study the effect of Lorentz contraction on the EOS properties is to 
incorporate them into transport models. The first steps in this direction have been made 
already in [26], but the approximation used in [26] is too crude. 

It might be more realistic to incorporate the developed approach into effective models of 
nuclear /hadronic matter [6, 7, 8, 9] and check the obtained EOS on a huge amount of data 
collected by the nuclear physics of intermediate energies. Since the suggested relativization 
of the VdW EOS makes it softer at high densities, one can hope to improve the description 
of the nuclear/hadronic matter properties (compressibility constant, elliptic flow, effective 
nucleon masses e.t.c.) at low temperatures and high baryonic densities [27]. 

Also it is possible that the momentum spectra of this type can help to extend the hydro- 
dynamic description into the region of large transversal momenta of hadrons (pr > 1.5 — 2 
GeV) which are usually thought to be too large to follow the hydrodynamic regime [28]. 

Another possibility to validate the suggested model is to study angular correlations of 
the hard core particles emitted from the neighboring regions and/or the enhancement of 
the particle yield of those hadrons occurring due to coalescence of the constituents with the 
short range repulsion. As shown above (also see Fig. 2), the present model predicts that the 
probability to find the neighboring particles with collinear velocities is higher than the one 
with non-collinear velocities. Due to this reason, the coalescence of particles with the parallel 
velocities should be enhanced. This effect amplifies if pressure is high and if particles are 
relativistic in the local rest frame. Therefore, it would be interesting to study the coalescence 
of any relativistic constituents with hard core repulsion (quarks or hadrons) at high pressures 
in a spirit of the recombination model of Ref. [29] and extend its results to lower transversal 
momenta of light hadrons. Perhaps, the inclusion of such an effect into consideration may 
essentially improve not only our understanding of the quark coalescence process, but also the 
formation of deuterons and other nuclear fragments in relativistic nuclear collisions. This 
subject is, however, outside the scope of the present work. 

As a typical VdW EOS, the present model should be valid for the low particle densities. 
Moreover, our analysis of the limit |i/T» 1 for fixed T leads to a surprisingly clear kinetic 
expression for the system's pressure (32). Therefore, it is possible that this low density result 
may provide a correct hint to study the relativistic analog of the dense packing problem. 
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Thus, it would be interesting to verify, whether the above approach remains valid for rela- 
tivistic quantum treatment because there are several unsolved problems for the systems of 
relativistic bosons and/or fermions which, on one hand, are related to the problems discussed 
here and, on the other hand, may potentially be important for relativistic nuclear collisions 
and for nuclear astrophysics. 
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In order to study the high pressure limit, it is necessary to estimate the excluded volume 
of two ellipsoids, obtained by the Lorentz contraction of the spheres. In general, this is quite 
an involved problem. Fortunately, our analysis requires only the ultrarelativistic limit when 
the mean energy per particle is high compared to the mass of the particle. The problem can 
be simplified further since it is sufficient to find an analytical expression for the relativistic 
excluded volume with the collinear particle velocities because the configurations with the 
noncollinear velocities have larger excluded volume and, hence, are suppressed. Therefore, 
one can safely consider the excluded volume produced by two contracted cylinders (disks) 
having the same proper volumes as the ellipsoids. For this purpose the cylinder's height in 
the local rest frame is fixed to be | of a sphere radius. 

Let us introduce the different radii R\ and R2 for the cylinders, and consider for the 
moment a zero height for the second cylinder hi = and non-zero height hi for the first 
one. Suppose that the center of the coordinate system coincides with the geometrical center 
of the first cylinder and the OZ-axis is perpendicular to the cylinder's base. Then the angle 
0„ between the particle velocities is also the angle between the bases of two cylinders. To 
simplify the expression for the pressure, the Lorentz frame is chosen to be the rest frame of 
the whole system. 

In order to estimate the excluded volume we fix the particle velocities and transfer the 
second cylinder around the first cylinder while keeping the angle 0„ fixed. The desired 
excluded volume is obtained as the volume occupied by the center of the second cylinder 
under these transformations. Considering the projection on the XOY plane (see Fig. 3. a), 
one should transfer the ellipse with the semiaxes R x = R2 cos (©„) and R y = R 2 around 
the circle of radius Ri. We approximate it by the circle of the averaged radius (Rxoy) = 
Ri + (R x + R y )/2 — Ri + i?2(l + cos (0„))/2. Then the first contribution to the excluded 
volume is the volume of the cylinder of the radius (Rxoy) and the height hi = CC\ of the 
cylinder OABC in Figs. 3. a and 3.b, i.e., 



Appendix A: Relativistic Excluded Volume 




(44) 
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Projecting the picture onto the XOZ plane as it is shown in Fig. 3.b, one finds that the 
translations of a zero width disk over the upper and lower bases of the first cylinder (the 
distance between the center of the disk and the base CA is, evidently, CD\ = R 2 \ sin (©„) |) 
generate the second conrtibution to the excluded volume 

vnih) = txR\ 2 R 2 \ sin (Q v ) | . (45) 

The third contribution follows from the translation of the disk from the cylinder's base 
to the cylinder's side as it is shown for the YOZ plane in Fig. 3.c. The area BB\F is the 
part of the ellipse segment whose magnitude depends on the x coordinate. However, one can 
approximate it as the quarter of the disk area projected onto the YOZ plane and can get a 
simple answer nR 2 \ sin (6„) |/4. Since there are four of such transformations, and they apply 
for all x coordinates of the first cylinder (the length is 2i?i), then the third contribution is 

vmih) = tiR\ 2 Ri\ sin (e„) | . (46) 

Collecting all the contributions, one obtains an estimate for the excluded volume of a 
cylinder and a disk 

v 2c (h 1 ) = n (#i + R 2 cos 2 (y)) h 1 + 2irR 1 R 2 (R 1 + R 2 )\sm(e v )\ . (47) 

The above equation, evidently, gives an exact result for a zero angle and arbitrary height of 
the first cylinder. Comparing it with the exact answer for 6„ = | 

vg (h 1: e v = I) = i?i (ttR 1 + AR 2 ) h + 2nR 1 R 2 (R 1 + R 2 ) , (48) 

one finds that the dominant terms (the second terms in (47) and (48)) again are exact, 
whereas the corresponding corrections, which are proportional to h\ , are related to each 
other as ~ 0.9897 (ratio of the approximate to exact values at R 2 = R±). Therefore, Eq. 
(47) also gives a good approximation for the intermediate angles and small heights. 

In order to get an expression for the non-zero height of the second cylinder we note 
that the expression for the excluded volume should be symmetric under the permutation 
of indexes 1 and 2. The lowest order correction in powers of the height comes from the 
contribution Vi(hi). Adding the symmertic contribution vj{h 2 ) to t>2 C (^i) (47), one obtains 
the following result 

v 2c (h 1 ,h 2 ) = 7i ( R l + R 2 cos 2 ( ^-\) h 1 + 7i(R 2 + R 1 cos 2 (^-)) h 2 



2 J J V V 2 

+ 2 77^2(^1 + ^2) sin . (49) 

The above expression gives an exact result for a zero angle and arbitrary heights of cylinders. 
It also gives nearly exact answer for Q v — | in either limit hi or h 2 — > 0. 

Choosing the heights to reproduce the proper volume for each of the Lorentz contracted 
spheres, one gets an approximation for the excluded volume of contracted spheres in ultra- 
relativistic limit 

v"r l (R u R 2 ) = ^(^ + ^cos 2 (^)) 2 + ^(i? 2 + i? 1 cos 2 (^)) 2 

+ 2nR 1 R 2 (R 1 + R 2 )sm(e v ) . (50) 
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C) 

Fig. 3. Relativistic excluded volume derivation for relativistic cylinder OABC and ultra- 
relativistic cylinder (disk) DC with radii Ri and R 2 , respectively. 6 is the angle between 
their velocities. Pictures a - c show the projections onto different planes. The transfer of the 
cylinder DC around the side of the cylinder OABC is depicted in Fig. 3. a. The solid curve 
DEF corresponds to the exact result, whereas the dashed curve corresponds to the average 
radius approximation (R X oy) = OA + (DC + BF)/2 = R ± + R 2 (l + cos (6))/2. 
The transfer of the cylinder DC = DC\ = AD 2 along the upper base of the cylinder 
OABC = ACC\A\ is shown in panel b. Its contribution to the excluded volume is a volume 
of the cylinder with the base AC = 2R\ and the height CD\ sin (0) = R 2 sin (©). A similar 
contribution corresponds to the disk transfer along the lower base of the cylinder A\C\. 
The third contribution to the relativistic excluded volume arises from the transformation 
of the cylinder DC = BBi = FB 1 from the upper base of the cylinder OABC = AB x O 
to its side, and it is schematically shown in Fig. 3.c. The area BB X F it/AR^ 1 sin (6) is 
approximated as the one quarter of the area of the ellipse BB\. 
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Fig. 4. Comparison of the relativistic excluded volume obtained by the approximative 
ultrarelativistic formula with the exact results. The left panel shows the quality of the 



approximation Vexcl 



,Urel 
'12 



(R, R) (50) to describe the excluded volume of two nonrela- 



tivistic spheres V2SP of the same radius R as a function of the spherical angle G. The right 
panel depicts the approximation to the excluded volume of the nonrelativistic sphere and 
disk. In both panels the solid curve corresponds to the exact result and the long dashed one 
corresponds to the ultrarelativistic approximation by two cylinders. The averaged ultrarel- 
ativistic excluded volume in the left panel is ^™g^ e ~ 1-065. The corresponding averaged 

value for the right panel is ^ VE y 2 ^ @ ~ 0.655, which should be compared with the exact value 



(Vexcl)® 
V2SP 



0.607. 



The corresponding 7 ? -factors (7^ = E(k q )/m q , q = {1, 2}) are defined in the local rest frame 
of the whole system for particles of mass m q . The last result is valid for < Q v < |, to use 



it for I < 0„ < it one has to replace 0^ 



7T 



e„ in (50). 



It is necessary to stress that the above formula gives a surprisingly good approximation 
even in nonrelativistic limit for the excluded volume of two spheres. For R2 = Ri = R 
one finds that the maximal excluded volume corresponds to the angle Q v = j and its 
value is max{t> ± 2 rel (R, R)} m y7ri? 3 , whereas the exact result for nonrelativistic spheres is 
t> 2s = y7ri? 3 , i.e., the ultrarelativistic formula (50) describes a nonrelativistic situation with 
the maximal deviation of about 10% (see the left panel in Fig. 4). 

Eq. (50) also describes the excluded volume v s d = 10 + 3?r R 3 for a nonrelativistic sphere 
and ultrarelativistic ellipsoid with the maximal deviation from the exact result of about 15% 
(see the right panel in Fig. 4). 

In order to improve the accuracy of (50) for nonrelativistic case, we introduce a factor 
a to normalize the integral of the excluded volume (50) over the whole solid angle to the 
volume of two spheres 



Nrel 



(i?l, i? 2 



^ „,Urel 

a v l2 



(Ri, R2 



a 



4tt (i?x + R 2 ) 



3jdQ v sin(e„) vYr\Ri,R2, 




(51) 



71=72=1 
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For the equal values of hard core radii and equal masses of particles the normalization 
factor reduces to the following value a ~ 1 Q g 54 , i.e., it compensates the most of the deviations 
discussed above. With such a correction the excluded volume (51) can be safely used for 
the nonrelativistic domain because in this case the VdW excluded volume effect is itself a 
correction to the ideal gas and, therefore, the remaining deviation from the exact result is 
of a higher order. 

It is useful to have the relativistic excluded volume expressed in terms of 3-momenta 

7i V 2 - R i l k ill k 2| / 7i V 

+ 2nR 1 R 2 (R 1 + R 2 ) 1 -^^ , (52) 

where Vo q denote the corresponding proper volumes Vo q = ^R q , Q — {1, 2}. 

For the practical calculations it is necessary to express the relativistic excluded volume 
in terms of the three 4- vectors - the two 4- momenta of particles, k qfJj , and the collective 
4- velocity w M = ^=f(l,v). For this purpose one should reexpress the gamma-factors and 
at least one of trigonometric functions in (50) in a covariant form 



|ki| 


N 


+ 1 


ki- 


k 2 | 


2R 2 


|ki 


|k 2 





7 ,= » +k ' A, C o 3( e„) = , ^K^-KW m 

m m \/((*!'«»,) 2 -m 2 )((*?i',.) 2 -m 2 ) 

Using Eq. (53), one can express any trigonometric function of 0„ in a covariant form. 
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